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ABSTRACT 

It is expected that specific globular clusters can contain up to a hundred of millisec- 
ond pulsars. These pulsars can accelerate leptons at the shock waves originated in 
collisions of the pulsar winds and/or inside the pulsar magnctospheres. Energetic lep- 
tons diffuse gradually through the globular cluster comptonizing stellar and microwave 
background radiation. We calculate the GeV-TeV 7-ray spectra for different models of 
injection of leptons and parameters of the globular clusters assuming reasonable, of the 
order of 1%, efficiency of energy conversion from the pulsar winds into the relativistic 
leptons. It is concluded that leptons accelerated in the globular cluster cores should 
produce well localized 7-ray sources which are concentric with these globular clusters. 
The results are shown for four specific globular clusters (47 Tuc, Ter 5, M13, and 
M15), in which significant population of millisecond pulsars have been already discov- 
ered. We argue that the best candidates, which might be potentially detected by the 
present Cherenkov telescopes and the planned satellite telescopes (AGILE, GLAST), 
are 47 Tuc on the southern hemisphere, and M13 on the northern hemisphere. We 
conclude that detection (or non-detection) of GeV-TeV 7-ray emission from GCs by 
these instruments put important constraints on the models of acceleration of leptons 
by millisecond pulsars. 

Key words: globular clusters: individual (47 Tuc, Ter 5, M15, M13) - pulsars: - 
radiation mechanisms: non-thermal - gamma-rays: theory 



1 INTRODUCTION 

Globular clusters contain ~10 5 -10 6 late type stars in a vol- 
ume with the half-mass radius of the order of a few parsecs. 
They are distributed spherically around our Galaxy at char- 
acteristic distances of ~10 kpc (e.g. Harris 1996). It is ex- 
pected that the typical globular cluster can contain up to a 
hundred of millisecond radio pulsars (MSP) which appear as 
a result of spinning up due to the accretion of matter from 
the low mass companion (the so called "recycling" model, 
e.g. Alpar et al. 1982). In fact, more than 100 radio pulsars 
with spin periods in the range of a few to a few tens of mil- 
liseconds, have been detected in 24 globular clusters. The 
largest observed samples are from Ter 5 (23 MSP) and 47 
Tuc (22 MSP) (e.g. Camilo & Rasio 2005). Many of these 
MSPs are inside the LMXB systems, but significant part is 
isolated probably due to the evaporation of the companion 
star. It is expected that MSPs produce relativistic magne- 
tized winds, which in the case of binary pulsars interact with 
the stellar winds. As a result, a shock wave is created on 
which particles can be accelerated to energies large enough 
for 7-ray production (e.g. Bigniami, Maraschi & Treves 1977, 
Kluzniak et al. 1988, Phinney et al. 1988). Only a part of 



E-mail: bednar@fizwe4.phys.uni.lodz.pl 



the pulsar wind is obscured due to the presence of the stel- 
lar companion. In most directions it propagates unconfined 
and interacts with the globular cluster medium or with other 
pulsar winds. Note, that in the case of very compact binaries 
pulsars can be permanently enshrouded by the evaporating 
material from the companion star (Tavani 1991). This class 
of "hidden" pulsars is difficult to detect and none of the 
objects of this type has been uniquely identified up to now. 

The 7-ray fluxes expected from specific binaries con- 
taining MSPs (e.g. PSR 1957+20, Arons & Tavani 1993) or 
from the population of single MSPs inside the globular clus- 
ters (Tavani 1993) are about the order of magnitude below 
the sensitivity of the EGRET telescope. Larger fluxes are ex- 
pected from the hidden MSPs in compact binaries (Tavani 
1991). In fact, only the upper limits on the 7-ray emission 
from the globular clusters are available, e.g. (1 — 2) x 10~ 7 ph. 
cm~ 2 s" 1 above 100 MeV (Michelson et al. 1994, Manand- 
har et al. 1996) and 8 x 10 -6 ph. cm -2 s _1 in energy range 
0.75 - 30 MeV (O'Flaherty et al. 1995). The above men- 
tioned upper limits already allow to put constraints on the 
number of such sources in globular clusters. More recently, 
7-ray emission from the inner MSP magnetospheres has been 
also estimated for specific GCs by Harding et al (2005). It 
has been concluded that 7-rays above 50 GeV should be ob- 
served from some specific individual MSPs and from some 
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GCs by the MAGIC, HESS, and VERITAS Cherenkov tele- 
scopes. 

The aim of this paper is to calculate the 7-ray emission 
in the GeV-TeV energy range expected from the whole pop- 
ulation of MSP in the globular clusters. However, in con- 
trast to previous models, we assume that 7-ray emission 
is produced by leptons which scatter the stellar and mi- 
crowave background radiation during their diffusion through 
the medium of the GC. Leptons are accelerated either on the 
shock waves produced in collisions of the MSP winds inside 
the GC or injected from the inner pulsar magnetopsheres. 
We show that leptons can be accelerated to high energies and 
are confined inside globular cluster for sufficiently long time 
for efficient up-scattering of the stellar and microwave back- 
ground radiation to GeV-TeV energy range. The available 
upper limits on the TeV 7-ray emission from the globular 
clusters, e.g. in the case of M 13 (Hall et al. 2003) and M 15 
(LeBohec et al. 2003), do not constrain strongly the possi- 
ble production of energetic particles in these type of objects. 
However, it is likely that GeV-TeV 7-rays can be detected 
from globular clusters by the present and future telescopes. 



2 PULSARS INSIDE GLOBULAR CLUSTERS 

The number of millisecond pulsars in the central part of 
typical globular cluster (GC) is expected to be of the order 
of a few tens up to a few hundred (e.g. Tavani 1993). In 
fact, 23 pulsars have been observed up to now in Ter 5 and 
22 pulsars in 47 Tuc. Most of these pulsars are localized 
inside the core of GCs. They have the average period of the 
order of a few milliseconds and surface magnetic fields of 
the order of ~10 9 G (see e.g. Camilo & Rasio 2005). The 
characteristic energy loss rate due to the pulsar wind from 
an isolated pulsar can be estimated by applying the formula 
for rotating magnetic dipole, 

L p = 1.2 x 10 35 B%P± A ergs" 1 , (1) 

where B p — 10 Bg G is the surface magnetic field of the 
pulsar, and P — 4P4 ms is the pulsar period. 

For typical radius of the core of GC, R c ~ 0.5 pc, the 
characteristic distances between pulsars in GC are R p w 
R c /Np^ 3 . For example, assuming the presence of 100 pul- 
sars in the cluster core, we obtain R p w 3.3 x 10 17 cm. Pul- 
sars produce relativistic strongly magnetized winds which, 
when colliding with other winds, create relativistic shocks. 
If the pulsar winds collide between themselves then typical 
distance from the pulsar to the shock can be estimated to 
be 

R sh ~ 0.5i? p w 0.5i? c /iVp 1/3 . (2) 

We estimate the magnetic field strength at such shocks by 
simple re-scaling from the pulsar (assuming i?~ 3 dependence 
in the inner pulsar magnetosphere and dependence in 
the pulsar wind), 

B sh « 3V^B p (R NS /RLc) 3 (RLc/R s h) 

« 2.5 x W~ 5 ^B 9 /Pi G, (3) 

where a is the magnetization parameter of the pulsar wind 
at the shock region, J?ns = 10 6 cm is the neutron star ra- 
dius, _Rlc = cP/2-k is the light cylinder radius. We consider 



the case of leptons injected by the pulsars and additionally 
accelerated in such relativistic shocks with the power-law 
spectra. The value of a has been estimated in the case of 
the wind from the Vela type pulsars at ~0.1 (e.g. Sefako & 
de Jager 2003) and from the Crab type pulsars at ~0.003 
(Kennel & Coroniti 1984). If the wind is terminated closer 
to the pulsar, then a may be of the order of ~1 (Contopou- 
los & Kazanas 2002). Applying this available range of a, we 
estimated the magnetic field at the pulsar wind shock region 
in the range B sh ~ 10~ 6 - 10" s G. 

The pulsar winds can also interact with the winds of 
the stars inside the GC. There are typically JV* ~2000 stars 
in the cores of the GCs within the spherical volume with 
the radius of R c = 0.5 pc. The average distance between 
these stars is then R* w R c /Nl /3 w 1.2 x 10 17 cm. Let 
us consider for simplicity solar type stars with typical mass 
loss rate of the order of M = 10~ 14 M per year and the 
wind velocity V* ~ 500 km s _1 . We determine the parame- 
ter n — Lp/cMVi,, which describes the pressure balance be- 
tween the pulsar and the stellar winds. For the above values 
it is equal to fx « 1.25 x 10 5 . Then, based on the simpli- 
fied model for the structure of the colliding winds (Girard 
& Wilson 1987), we estimate the radius of the volume dom- 
inated by the stellar wind on R^ = R*/(l + ^/Ji), i.e. equal 
to Rd « 2.8 x 10~ 3 i?* w 3.4 x 10 14 cm for a solar type star. 
Therefore, we conclude that stellar winds are confined in a 
small region around the stars and the pulsar winds collide 
mainly between themselves. 

Globular clusters can contain typically a few hundred 
thousand stars with characteristic masses of the order of 
~ 1 Mq in the volume of a few parsecs. Since their mass to 
luminosity ratio is close to ~ 2 (see e.g. recent analysis for 
M15 by van den Bosch et al. 2006), these stars create very 
strong background radiation field which can act as a target 
for relativistic leptons injected by the millisecond pulsars. 
Based on the observed luminosity of specific globular cluster, 
Lgc, and on the density profile for the distribution of the 
stars inside the cluster (based on the Michie-King model 
(Michie 1963), see also Kuranov & Postnov 2006), 

C 1, R<R C 
D{R) = <^ (Rc/R) 2 , R c < R < Rh (4) 
[ (R c Rh) 2 /R 4 , Rh < R < Rt, 

we calculate the energy density of stellar photons inside the 
cluster core from the formula, 

_ L GC 18fl? - 3^/6RfR c - 2R\ 

Urad - cr* 6{ vmm-2RZ) ' ( ) 

where R c is the GC core radius; Rt is its tidal radius, 
and the half mass radius, Rh, is related to those both by 
Rh = y // 2R c Rt/3. For typical parameters, R c = 0.5 pc, 
Rt = 50 pc, and L G c = 10 5 L Q = 3 x 10 3S erg s _1 , we ob- 
tain C/ ra d ~ 300 eV cm~ 3 . The energy density of this stellar 
radiation field clearly dominates over the energy density of 
the microwave background radiation (MBR), Umbr. = 0.25 
eV cm -3 . It can also dominate over the energy density of 
the magnetic field inside the GC provided that the mag- 
netic field is not stronger than Bgc ~ 10~ 4 G. This limit 
has been obtained assuming the equipartition between the 
energy density of stellar photons and the energy density of 
the magnetic field. Since such strong magnetic field inside 
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the GCs seems unlikely (see estimate given by Eq. [3] for the 
region of the pulsar wind shock), we conclude that leptons 
injected by the pulsars in the GC cores lose energy mainly 
on IC scattering of the stellar and microwave radiation. 



3 ACCELERATION OF LEPTONS INSIDE 
GLOBULAR CLUSTERS 

Let us estimate the maximum energies which leptons (e ) 
can reach in such pulsar wind shocks. The obvious limit is 
introduced by the condition that the Larmor radius of ac- 
celerated particles can not be greater than the dimension of 
the shocks, i.e. Rl = cp/(eB) < R s h, where c is the velocity 
of light, e is the charge of an electron, p is the momentum 
of the particle, and B is the magnetic field strength at the 
shock. For the above estimated magnetic field at the shock, 
we get 

S S T X < 1.4 x W 3 ^B g /Pi TeV. (6) 

Additional limits can be constrained by the energy losses of 
leptons on the synchrotron and inverse Compton (IC) pro- 
cesses during their acceleration. In order to estimate these 
limits we compare the acceleration rates of leptons with their 
energy loss rates on the mentioned processes. The accelera- 
tion rate is independent on the energy of lepton E provided 
that so called acceleration coefficient £ is constant. It can be 
written as, 

= CcS/^l ~ lO 13 ^ eVs" 1 , (7) 

\ OA* J acc 

and the corresponding acceleration time scale, 

tacc = E/(dE/dt) acc » 10 5 £ Tc v/(eB-6) s. (8) 

where E = 1-BtoV TeV is the energy of leptons, B s h = 
10 -6 B-6 G is the magnetic field strength at the shock, and 
c is the velocity of light. 

On the other hand, the synchrotron energy loss rate is, 

= ^^Tf/B7c « 8.3 x 10- 4 B s 2 h7c 2 eVs" 1 , (9) 

where 17b is the energy density of magnetic field, <tt is the 
Thomson cross section, and 7 C is the Lorentz factor of lep- 
tons. By comparing Eqs. [7] with [9] and applying Eq.|3]we get 
the limit on the maximum energy of accelerated leptons due 
to the synchrotron losses, 

E™T « 4 x 10 2 P 4 d / 2 2 /(^ 1/4 -B9 /2 ) TeV. (10) 

where £ = 0.01£_2- This limit is more restrictive than the 
limit due to the dimension of the shock (-E"yn X < E™*) 
provided that the following condition is fulfilled, 

P 4 < l.Ba^B^/eJ?- (11) 

Leptons accelerated at the shock are also immersed in 
the strong radiation field coming from the stellar population 
within the GC. Above we estimated the average density of 
stellar photons in the core of GC on C/ ra d ~ 300 eV cm -3 . 
This allows us to estimate the energy loss rate on ICS in the 
Thomson (T) regime, 

(?)^ = ^a T C/ rad7 2 ^7.6xl0- 12 7c 2 eVs" 1 . (12) 
The Thomson regime is valid only for leptons with energies 



below Et/kn ~ 100 GeV. By comparing Eqs.[7lwith ll2l and 
applying the limit -Et/kn, we obtain that the saturation of 
acceleration of leptons can occur due to IC losses in the 
Thomson regime provided that £ < £ C rit = 10 -7 /B-6- How- 
ever, this case is not specially interesting from the point of 
view of TeV 7-ray production. Note that the acceleration of 
particles at relativistic shocks (defined by £) is usually con- 
sidered as much more efficient. For example, observations of 
the TeV 7-ray emission from the Crab Nebula show that the 
acceleration of particles occurs during a few Larmor cycles. 
Thus, £ has to be of the order of ~ 0.1 — 1. Therefore, here 
we consider the case of saturation of lepton acceleration in 
the Klein-Nishina (KN) regime. The energy loss rate in the 
KN regime depends only logarithmically on energy of lep- 
tons, i.e. depends very weakly on energy of leptons. So then, 
if the acceleration coefficient is significantly larger than £ cr it 
(as considered in this paper for the pulsar wind relativis- 
tic shocks) the saturation of acceleration of leptons might 
be caused either by synchrotron energy losses or dimension 
of the shock (depending on the parameters of the specific 
pulsar, see Eq. Hip . 

However, the limit on the maximum energies of accel- 
erated leptons can be also given by the advection time scale 
of leptons along the surface of the shock. The velocity of 
the plasma flow downstream of the relativistic shock can be 
estimated at v a dv ~ c/3. Then, the advection time scale is, 

tadv = iW«adv = 0.5 J R c /(Af p 1/3 Uadv). (13) 

where is defined by Eq.0 By comparing the acceleration 
time scale of leptons (Eq. [8} with their advection time scale 
along the shock (Eq. 113ft . we obtain another limit on the 
maximum energies of leptons, 

E™<?v = 1.5^_2-B-6-Rpc/(A r i / 3 i | c) TeV, (14) 

where R c = li? pc pc, jV p = lOOiVioo, and v c = u a dv/c. 
For likely parameters of the considered scenario, _B_e = 1, 
R pc = 0.5, TVioo = 0.3, and v c — 0.3, we estimate the max- 
imum energies of accelerated leptons due to their advec- 
tion along the shock at E£&? « 370£ TeV. For reasonable 
acceleration efficiency at relativistic shock of the order of 
£ ~ 0.1 — 0.01, this maximum energies of leptons are limited 
to the range E™* ~ 4 — 40 TeV. Therefore, we conclude 
that the acceleration process of leptons is limited by their 
advection process along the surface of the shock. 

Based on the above general considerations, we conclude 
that the maximum energies of leptons accelerated at the 
shocks are likely to be limited by advection process of lep- 
tons along the shock and not by the dimension of the shock 
in the pulsar wind -E^ ax or the energy losses on the syn- 
chrotron or IC scattering processes. Leptons can be acceler- 
ated inside GC to energies which allow them to comptonize 
the stellar and MBR photons to the TeV 7-ray energy range. 

Leptons injected with different energies in the central 
core of the globular cluster diffuse gradually in the outward 
direction. We consider the central region with the radius R^ 
in which half of the mass of the GC is contained. Assuming 
the Bohm diffusion, the average time spend by leptons with 
energy E inside this region is, 

tdiff = Rl/D diS , (15) 

where Ddiff = Rlc/3 is the diffusion coefficient, _Rl = 
ep/eBcc « 3 x 10 15 £ T cv/B- 6 cm is the Larmor radius of 
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Figure 1. The density of stars (figure a), and density and energy 
density of stellar photons (b) as a function of distance R from 
the center of a typical globular cluster with the mass 10 Mq, 
the core radius R c = 0.5 pc, the half mass radius i?h = 4 pc, 
the tidal radius Rt = 50 pc. The stellar density profile is defined 
by Eq. [4] The photon densities and energy densities of stellar 
photons are shown by the thick solid and dashed curves and the 
corresponding values for the microwave background radiation are 
shown by the thin horizontal lines. 

leptons in the magnetic field of the GC, and Bqc = 10 _6 i?_6 
G is the average magnetic field in the GC. 

In order to answer the question whether leptons can lose 
significant part of their energy during the diffusion through 
the thermal radiation field inside the GC we estimate the col- 
lision rate on IC process in the Thomson regime for supposed 
parameters inside the GC. Let us use the following parame- 
ters, EtcV = 0.1, Bgc — 3 x 10 -6 G, i?h = 3 pc, and the av- 
erage density of thermal photons inside the GC, N p h = 100 
cm -3 . Then, i? L = 10 14 cm, D diff = 10 24 cm 2 s _1 , and 
idift = 8x 10 13 s. For this diffusion time scale, the collision 
rate of leptons is of the order of AT co i = ctdiacnNph ~ 160, 
i.e. leptons might interact up to a hundred times before leav- 
ing the globular cluster. Note that the diffusion time and the 
ICS cross section in the KN regime are inversely propor- 
tional to energy of leptons. Therefore, leptons with energies 
up to ~ 1 TeV should interact frequently with stellar pho- 
tons producing 7-rays with comparable energies. However, 
leptons with energies above ~ 1 TeV have chance to escape 
from the radiation field of the GC without interaction with 
stellar photons. On the other hand, these high energy lep- 
tons interact frequently with the MBR producing 7-rays in 
the energy range from tens MeV up to hundreds GeV. In 
the next section, we perform more detailed calculations of 
the 7-ray spectra produced in the process considered above 
by assuming injection of leptons with the mono-energetic 
and the power-law spectra with different spectral indexes, 
injection efficiencies, and high energy cut-offs. 



4 GAMMA-RAY PRODUCTION INSIDE 
GLOBULAR CLUSTERS 

Let us assume that the part, 77, of the wind energy of pulsars 
inside the GC is converted into relativistic leptons (so called 
energy conversion efficiency). Then using Eq. [1] the power 
injected in relativistic leptons is, 

L c = rjN p L p = 1.2 x 10 35 r/A/pB|P 4 - 4 erg s _1 . (16) 



The spectrum of leptons accelerated at the shocks can be 
approximated by a single power-law with the spectral index 
a between -E m i n and -Emax- These leptons produce 7-rays 
as a result of ICS of thermal and MB photons. In order 
to perform numerical calculations of the 7-ray spectra, first 
we have to define the radiation field inside GC as a func- 
tion of the distance from its center. The density profile for 
the distribution of stars inside GC (Eq. |4j is normalized 
to the total number of stars inside GC which is typically 
in the range AT* ot = 10 5 - 10 6 Mq. The density of stars 
as a function of the distance is N ir (R) = AD(R)dR, where 
A = Nl ot /[2R c - 2R 2 c /(3R h ) - R*Rl/(3R%)] « Nl ot /(2R c ) 
and D(R) is given by Eq. [4] As a first approximation, it 
is assumed that the average star inside GC has the basic 
parameters close to that of our Sun, i.e. the surface tem- 
perature 6000 K and the radius 7 x 10 10 cm (note that the 
mass to light ratio for typical GC is close to ~2). The den- 
sity of stellar photons at a specific distance from the center 
of GC is calculated by integration over the whole distribu- 
tion of stars inside the GC. The examples of the density 
profiles for distribution of stars and the density and energy 
density of stellar photons inside the GC with the luminos- 
ity Lgc = 10 5 Mq are shown in Fig. [1] Inside the globular 
cluster core, defined by R c = 0.4 pc, the density of photons 
drops slowly with the distance from the core reaching the 
value of approximately a factor of two lower at R c than at 
the center of the core. Outside the half mass radius, defined 
by Rh — 4 pc, density of photons drops almost proportion- 
ally to oc R~ 2 . For comparison we also show the photon and 
photon energy densities of MBR. Note, that in the case of 
considered GC, the photon density of the MBR dominates 
over density of stellar photons. In contrast, the photon en- 
ergy density of MBR is much lower than the stellar photon 
energy density in the whole volume of the GC. Therefore low 
energy leptons lose mainly energy on scattering stellar ra- 
diation in the Thomson regime. However leptons with large 
energies lose energy mainly by scattering the MBR in the 
Thomson regime since their scattering of stellar radiation is 
strongly suppressed due to the KN effect. 

4.1 Injection of mono-energetic leptons 

We calculate the 7-ray spectra produced by leptons which IC 
up-scatter the stellar photons and the MBR. For considered 
energies of leptons the scattering process of stellar radiation 
can occur in the T or the KN regimes. We neglect the syn- 
chrotron energy losses of leptons for the considered range 
of their energies and values of the magnetic field inside the 
GCs. It is assumed that leptons are injected homogeneously 
inside the core of GC. Since diffusion of leptons from the core 
occurs in the variable thermal radiation field, the impor- 
tance of the scattering process in both considered radiation 
fields (thermal photons and MBR) can differ significantly. In 
order to calculate the 7-ray spectra produced by leptons dif- 
fusing outward the center of GC we apply the Monte Carlo 
method. The optical depth for leptons on the ICS is deter- 
mined during their gradual diffusion in variable radiation 
field. Therefore, the distance from the core of GC at which 
lepton with specific energy E produces 7-ray photon is sim- 
ulated numerically. For known location of the interaction 
place inside the GC, we simulate the energy of 7-ray photon 
produced in the IC process. The energy of primary lepton 
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Figure 2. Differential 7-ray spectra (multiplied by the energy 
squared) produced in the IC scattering of the stellar and MBR 
radiation by mono-energetic leptons with energies E = 10 GeV 
(triple-dot-dashed curve), 10 2 GeV (dotted), 10 3 GeV (dot- 
dashed), 10 4 GeV (dashed), and 10 5 GeV (solid), after normal- 
ization to a single lepton per second. Leptons are injected by the 
millisecond pulsars in the core of GC and diffuse in the outward 
direction in the magnetic and radiation field created by the stars 
and the MBR. The magnetic field strength is assumed equal to 
Bqq = 10 -6 G. The diffusion process of leptons is followed up to 
100 pc from the core of GC. The total stellar luminosity of GC is 
equal to L GC = 10 5 L Q (a) and 10 6 L o (b). 



is reduced by the energy of produced 7-ray and the process 
of simulation is repeated for the next interaction place and 
energy of produced 7-ray photon. The propagation process 
of leptons is considered up to the distance of 100 pc from 
the center of the GC. At this distance, the contribution of 
stellar photons becomes unimportant. All secondary 7-rays 
produced by leptons are sorted in specific energy ranges in 
order to obtain the spectrum of 7-rays produced inside the 
GC. The method described above allows us to calculate the 
spectra of produced 7-rays as a function of distance from 
the center of GC. Therefore, we can estimate the extent of 
the specific 7-ray source on the sky related to the GC. 

In order to perform specific calculations we assume that 
leptons are injected homogeneously in the core of GC. Dif- 
ferential 7-ray spectra (multiplied by the energy squared), 
calculated for mono- energetic leptons injected at the core of 
GC, are shown in Fig. [2] for the total luminosity of the GC 
Lgc = lO 5 I/0 (figure a) and lO 6 Z/0 (b) after normalization 
to a single injected lepton per second. For leptons with low 
energies 10 2 GeV), the IC scattering of stellar photons 
and MBR occurs only in the T regime. As we noted above 
the photon density of MBR dominates over the density of 
stellar photons for Lgc = W 5 Lq, but the energy density of 
the MBR is significantly lower. Therefore, such leptons lose 



energy mainly by scattering stellar photons. For the total 
luminosity of GC equal to Lgc = 10 6 Lq both the photon 
density and energy density of stellar photons dominate over 
the MBR photon and energy density. However, when the en- 
ergies of injected leptons are larger, the situation becomes 
more complicated due to the influence of the ICS in the KN 
regime. Then, for sufficiently large energies of leptons, the 
scattering process of MBR can dominate the energy losses 
of leptons. 

Let us at first consider the case of the GC with lumi- 
nosity Lgc = 1O 5 L0. It is clear from our calculations that 
scattering of MBR can dominate energy loss process of lep- 
tons. Already for leptons with energies 10 3 GeV, the power 
which goes into the 7-rays due to scattering of MBR in the 
T regime and the stellar radiation in the KN regime be- 
comes comparable (see the 7-ray spectrum marked by the 
dot-dashed curve in Fig. [2^, with the characteristic lower en- 
ergy 'bump' due to IC scattering of MBR in the T regime 
and the characteristic 'peak' at the highest energies due to 
the IC scattering of the stellar radiation in the KN regime) . 
For the largest considered energies of leptons (10 5 GeV) only 
the 7-ray spectrum produced in the IC scattering process of 
the MBR is clearly seen but the scattering of stellar photons 
is completely suppressed due to the KN regime. However, 
when the density of stellar photons dominates over the den- 
sity of the MBR (the case of the GC with the total luminos- 
ity equal to Lgc = 1O 6 L0), then contribution to the 7-ray 
spectrum from leptons which scatter also stellar photons in 
the KN regime becomes evident up to the largest considered 
energies of leptons (see characteristic bumps at the end of 
7-ray spectrum due to scattering in the KN regime, Fig.[2p). 
Note, that in case of such luminous GCs, the 7-ray spectra 
produced by leptons with energies above 10 4 GeV are clearly 
dominated by photons produced in the scattering process of 
the MBR. For leptons with energies 10 3 GeV, the contri- 
bution from scattering of the stellar and the MBR photons 
becomes comparable but for leptons with energies below 10 2 
GeV, 7-ray spectrum is formed mainly by scattering the stel- 
lar radiation since density of these photons in the core of the 
GC dominates over the density of the MBR by a factor of 
~3. 

4.2 Injection of leptons with the power 
law-spectra 

As a next step, we analyze the features of the 7-ray spectra 
produced by leptons with the power-law spectrum. As we 
discussed above (Sect. 3), leptons are accelerated with the 
power-law spectrum at the relativistic shocks which are cre- 
ated in collisions of the pulsar winds inside the core of GCs. 
The maximum energies of leptons, B m „, are determined by 
the radiation and escape processes of leptons and by the ef- 
ficiency of their acceleration. For illustration purposes, we 
consider in this section the cases with L max = 3, and 30 TeV 
which are consistent with the above derivations (Sect. 3). 
Since details of the acceleration process are not precisely 
known, we consider a range of the injection spectra with 
the spectral indexes between a=2-3. On the other hand, we 
also fix the low energy cut-off in these power-law spectra on 
-Emm = 1 GeV or 100 GeV. This low energy cut-off might 
correspond to energies of leptons injected into the pulsar 
wind from the inner magnetopsheres of the millisecond pul- 
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sars. The models of the radiation processes in the inner mag- 
netospheres of MSPs show that if the acceleration of leptons 
occurs as expected in terms of a space charge limited flow 
model (developed by Muslinov & Tsygan 1992), then their 
acceleration can be limited by the curvature energy losses to 
energies ~1-10 TeV (e.g. Bulik, Rudak & Dysk 2000; Luo, 
Shibata, Melrose 2000; Harding, Usov & Muslimov 2005). 
In fact, this possibility is investigated in this paper in the 
case of injection of mono-energetic leptons (see Sect 4.1 and 
5.5). However, the processes occurring in the inner magne- 
tospheres of MSPs are not well known. If the temperature 
of the polar cap region is sufficiently large, above ~ 3 x 10 6 
K, then plenty of e pairs can appear in the acceleration 
region due to 7-ray production in ICS process and 7-ray 
photon absorption in the magnetic field or in two photon 
collisions (Harding & Muslinov 2002; Harding, Muslinov & 
Zhang 2002). These e pairs may saturate the electric field 
in the acceleration region. As a result, leptons may escape 
from the inner magnetosphere with relatively low energies. 

The power-law spectrum of leptons is normalized to the 
part of power of the millisecond pulsars winds present inside 
the core of GC, L e (see Eg . 116(1 . assuming the average pulsar 
parameters, Bg = 1 and P4 = 1. The energy conversion ef- 
ficiency of the pulsar energy loss rate to relativistic leptons 
is taken equal to 77 = 1%, and the number of millisecond 
pulsars in the GC core is equal to JV P = 30. In order to per- 
form the Monte Carlo simulations, we divide the spectrum 
of leptons into the narrow energy bins (E, E + AE) since 
the diffusion and ICS processes of leptons strongly depends 
on their energy. 

To understand better the 7-ray production process, we 
calculate the differential 7-ray spectra produced by leptons 
in the IC scattering process of the MBR only and stellar ra- 
diation only (see dotted and dashed curves in Figs. [3^, and 
b). We also show the 7-ray spectra produced in the pres- 
ence of both radiation fields (solid curves in Figs. [3J. The 
power-law injection spectrum of leptons with spectral in- 
dex equal to a = 2.1, extending between i? m m = 100 GeV 
and i5 max = 3 x 10 4 GeV, is investigated. We follow diffu- 
sion of leptons up to the distance of 100 pc from the core 
of GC. Note that since the scattering of MBR by leptons 
occurs only in the T regime, the 7-ray spectrum, obtained 
in the case of the MBR only, has the shape characteris- 
tic for completely cooled leptons and the power-law 7-ray 
spectrum with spectral index close to (a + 2)/2 = 2.05. 
The 7-ray spectra obtained from IC scattering of stellar ra- 
diation only are completely different. As already discussed 
in Sect. 4.1, leptons with energies above ~0.1 TeV inter- 
act with stellar photons in the KN regime. Some of them 
can escape from the region with radius 100 pc around the 
core of GC without interaction. Therefore, the highest en- 
ergy part of the 7-ray spectrum becomes steeper and this 
effect is stronger for less luminous GCs. The 7-ray spectra 
obtained in the case of both photon fields (stellar and MBR) 
are combination of the above described spectra. In the case 
of the GCs with low luminosity (Lgc = 10 5 L©, Fig. [3^,), 
leptons scattering the MBR contribute mainly to the lower 
energy part of the 7-ray spectrum (see differences between 
dashed and full curves in Fig. [3^0 . If the luminosity of GC 
is large (Lgc = 10 6 L Q , Fig. [SJj), the contribution to the 
7-ray spectrum at low energies from scattering of the MBR 
is negligible since the energy and photon density of stellar 
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Figure 3. Differential 7-ray spectra (multiplied by the energy 
squared) produced by leptons with the power law spectrum and 
spectral index ct = 2.1 (extending between E mln = 10 2 GeV 
and i?m ax — 3 x 10 4 GeV) in the presence of only: MBR (dot- 
ted curve), stellar radiation (dashed), and both MBR and stellar 
radiations (solid). The spectrum of leptons is normalized to the 
power of 30 millisecond pulsars (see Eq, ll6l l. the energy conversion 
efficiency is assumed equal to tj = 0.01, and the average pulsar 
parameters are Bg = 1 and P4 = 1. The results for the the total 
stellar luminosity of the GC equal to Lqq = 10 5 Lq are shown in 
figure (a) and for 10 6 Lq in figure (b). The magnetic field strength 
inside the GC is fixed on B GC = 10~ 6 G. 



radiation completely dominates over the energy and stellar 
density of the MBR. The scattering of the MBR by leptons 
with the highest energies has some small effect on the pro- 
duced spectrum of 7-rays (shift to lower energies in respect 
to the 7-ray spectra from scattering of stellar radiation only 
due to the additional energy losses on the MBR). 

We also consider the dependence of the 7-ray spectra on 
the parameters describing the injection spectrum of leptons 
(Figs. HJi and b). For the GC with luminosity Lgc = 10 5 
M©, we calculate the 7-ray spectra for two spectral indices 
of leptons (a — 2.1 and 3), two low energy cut-offs (L m i n = 1 
and 100 GeV), and two high energy cut-offs (E m&K = 3 and 
30 TeV). As in the previous cases, the diffusion of leptons 
up to 100 pc from the center of GC is considered. The 7- 
ray spectra peak at the GeV or the TeV energy ranges de- 
pending on -E m i n . Therefore, the broad band 7-ray observa- 
tions, through GeV- TeV energy range, with the GLAST and 
Cherenkov telescopes (MAGIC, HESS, VERITAS) should 
confirm (or disprove) models considered here for different 
spectra of injected leptons. The detection (or non-detection) 
of 7-ray emission from GCs should put a limit on the pa- 
rameters of the acceleration models of leptons in inner mag- 
netopsheres of millisecond pulsars and/or at the pulsar wind 
relativistic shocks. 
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Figure 4. The differential spectra of 7-rays (multiplied by the 
energy squared) produced by leptons with the power law spec- 
trum and spectral indices a = 2.1 (solid curves), and 3 (dashed). 
The spectrum of leptons extends between -E m in = 100 GeV (thick 
curves) or 1 GeV (thin curves) and £ max = 30 TeV (a) or 3 TeV 
(b). The spectra of leptons are normalized as in Fig. [3] The GC 
with the luminosity Lgc = 10 Mq is considered. 



4.3 Dependence on the distance from the center 
of Globular Cluster 

Although globular clusters are typically at large distances, 
the regions on the sky from which 7-rays can arrive to the 
observer at the Earth can be quite large, e.g. at the typical 
distance of 5 kpc and the half mass radius 5 pc the GC 
has angular dimension of ~ 3'. Therefore, it is important to 
find out whether the 7-rays produced by leptons diffusing 
outside the central region of GC can create a point like or 
an extended source for the Cherenkov telescopes. In order to 
check this, we perform calculations of the 7-ray spectra as 
a function of distance from the center of the GC in case of 
their production by leptons with specific injection spectrum. 
The 7-ray spectra produced within the central volume with 
the radius 1 pc, 3.3 pc, 10 pc, 33 pc , and 100 pc are shown 
in Fig. [5] They were calculated for two different values of the 
magnetic field inside the GC, -Bgc = 10 -5 G (a) and 10 -6 
G (b). It is clear that in the case of the magnetic field equal 
to 10 G, most of the 7-rays are already produced within 
~ 3 pc from the center of the GC. The spectrum produced 
within ~1 pc is only a factor of 2 lower than within ~3 
pc at energies above a few hundred GeV. In the case of 
the GC with the magnetic field equal to 10 -6 G, the 7-ray 
emission is more spreaded-out. However, still most of the 
emission (~50%) comes from the central region with the 
radius of ~ 3 pc. Therefore, we conclude that the extent of 
the 7-ray source, produced by leptons injected within the 
GCs, is relatively small provided that the magnetic fields 



Figure 5. The 7-ray spectra (multiplied by the energy squared) 
produced by leptons with the power law spectrum and spectral in- 
dex 2.1 (between 100 GeV to 30 TeV) within the central region of 
the GC with the radius R<1 pc (triple-dot curve), <3.3 pc (dot- 
ted), <10 pc (dot-dashed), <33 pc (dashed), and 100 pc (solid). 
The magnetic field inside the GC region is fixed on Bqq = 10 — 5 G 
(a) and 10" 6 G (b). The specific GC has the core radius R c =0.5 
pc and the half mass radius R h = 4 pc. The GC with the lumi- 
nosity Lgc = 10 5 Mq is considered. 

inside the GC cores are of the order of 10" 6 - 10~ 5 G. For 
typical distance to the GC (of the order of 5 kpc), the 7- 
ray source with dimension of 3-5 pc should have an angular 
size of ~2'-3', i.e. it should be a point-like source for the 
present Cherenkov telescopes that have the typical angular 
resolution of this order. Since the extent of the 7-ray source 
is not very sensitive to the considered value of the magnetic 
field strength inside the GC, the more complicated model for 
the magnetic field structure inside the GC (e.g. dependence 
of the magnetic field on the distance from the center of GC) 
should not give very different results. 



5 EXPECTED GAMMA-RAY FLUXES FROM 
SPECIFIC OBJECTS 

As an example, we perform calculations of the 7-ray fluxes 
expected on the Earth from 4 globular clusters (selected 
from the catalog of globular clusters, Harris 1996). Inside 
these clusters many millisecond pulsars have been already 
detected and a large number of millisecond pulsars is ex- 
pected (see Camilo & Rasio 2005; Tavani 1993): 47 Tuc (22 
MSP detected and ~ 260 MSP expected), Ter 5 (23, and 
~ 230 MSP), M15 (8 and ~ 117 MSP), and M13 (5 and 
~ 102 MSP). Moreover, these clusters have different declina- 
tions. M15 and M13 can be observed at small zenith angles 
by the northern sky telescopes (MAGIC, and VERITAS), 
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and Ter 5 and 47 Tuc can be observed by the southern tele- 
scopes (HESS and CANGAROO). M15 and Ter 5 can be 
observed simultaneously by the northern and southern sky 
telescopes at small or large zenith angles what allows to 
cover the 7-ray energy range between ~100 GeV and a few 
tens TeV. This is important if one wants to extract possible 
variable contribution to the 7-ray flux from e.g., X-ray bi- 
naries or intermediate mass black holes also expected inside 
the GCs. 

For every considered here cluster, we calculate the stel- 
lar radiation field as a function of distance applying their 
known parameters, i.e., dimension and total luminosity. In 
all simulations we consider the propagation of leptons in- 
side the inner region around the center of GC with the ra- 
dius f?max = 10 pc. As we have shown above, most of the 
7-rays are produced in a relatively small region, provided 
that the magnetic field strength inside the cluster is larger 
than B GC = 10" 6 G. The 7-ray spectra are calculated for 
all these 4 clusters by assuming that the power injected in 
relativistic leptons is equal to L c ± = 1.2 x 10 35 erg s _1 . This 
power might be supplied by e.g., N p — 100 MSPs with the 
characteristic periods of 4 ms and surface magnetic fields 
of 10 9 G, and assuming the efficiency of energy conversion 
from the pulsar winds into the relativistic leptons equal to 
77 = 0.01 (see Eq. [16}. Note that the periods of MSPs de- 
tected in 47 Tuc are limited to a relatively narrow range 
around the average value of ~4 ms (see Fig. 3 in Camilo 
& Rasio 2005). We assume that leptons are accelerated at 
the pulsar wind shocks with the power-law spectra (oc E~ a ) 
and spectral indices in the range a — 2.1 — 3.0, between 
E mia = 1 GeV or 100 GeV and £ max = 3 TeV or 30 TeV. 
The applied energy range for spectrum of leptons has been 
motivated in Sect. 4.2. The results of calculations for specific 
GCs mentioned above are discussed in a more detail below. 

5.1 47 Tuc (NGC 104) 

47 Tuc is one of the closest globular clusters, at the distance 
of only 4.5 kpc. It belongs to the class of the most mas- 
sive GCs in our Galaxy. It has a core radius 0.44 arc min 
(R c — 0.58 pc), a half mass radius 2.79 arc min (Rh = 3.67 
pc), and the total stellar luminosity 7.5 x 10 5 L©. 22 mil- 
lisecond pulsars (with the periods in the range 2-6 ms) have 
been detected up to now within this cluster (e.g. Camilo 
et al. 2000) . However, much larger number of MSPs, of the 
order of ~ 200 is predicted in 47 Tuc (see Tavani 1993, or 
more recently Ivanova, Fregeau & Rasio 2004). 

The 7-ray spectra produced by leptons diffusing 
through the GC are shown in Fig. (6^,b, together with the 
sensitivities of the present and future 7-ray telescopes in 
the GeV and TeV energy range. For applied parameters of 
the acceleration scenario (JV P = 100 MSP and r\ = 0.01), 
these 7-ray fluxes should be easily detected by the GLAST 
telescope provided that the spectrum of leptons extends 
below 100 GeV for both considered spectral indices. The 
Cherenkov telescopes can clearly detect 7-ray emission from 
this GC for the considered model parameters even if only 
~20 already observed MSPs are present inside 47 Tuc. The 
only exception is the case of a steep spectrum of leptons 
extending significantly below 100 GeV (see e.g., results for 
spectral index a — 3.0). 

In case of non-detection by the HESS telescope, the 



product N p ■ 77, which determines the flux level at 7-ray ener- 
gies, should be ^0.01 for leptons with the spectrum extend- 
ing between 100 GeV and 30 TeV, ^0.03 for the spectrum 
between 1 GeV and 30 TeV and spectral index equal to 2.1, 
and SC2 for the spectral index equal to 3. The upper limits 
on the product N p ■ rj, which can be obtained by the GLAST 
telescope, are usually not so restrictive (see details in Table 

The available upper limits on the 7-ray flux in the MeV 
and GeV energy ranges from the direction of this cluster 
(mentioned in the introduction) do not put strong con- 
straints on the number of the millisecond pulsars and effi- 
ciencies of lepton acceleration (see e.g. Michelson et al. 1994, 
O'Flaherty et al. 1995, Manandhar et al. 1996). For example, 
in order not to exceed the EGRET upper limit, the product 
N p • r\ can not be larger than few, i.e. still a few hundred 
MSP might be present inside 47 Tuc provided that energy 
conversion efficiency is of the order of 77 = 0.01. 



5.2 Ter 5 

Ter 5 has a core radius 0.18 arc min (R c = 0.54 pc) and 
a half mass radius 0.83 arc min (Rh = 2.5 pc) at the dis- 
tance of 10.4 kpc. Its total stellar luminosity is 1.5 x 10 
L0. It contains 23 millisecond pulsars, the largest number 
detected so far within a single GC (Camilo et al. 2000). Ter 
5 has been observed by the EGRET telescope, but derived 
upper limit does not constrain significantly the population 
of the millisecond pulsars and the acceleration efficiencies of 
leptons (Michelson et al. 1994) . The 7-ray spectra calculated 
for this GC and variety of considered models for lepton ac- 
celeration are shown in Fig. [7] Due to the smaller luminosity 
and larger distance, these 7-ray spectra are on a significantly 
lower level than those previously estimated for 47 Tuc. They 
are close to the detection limit of the GLAST and MAGIC 
telescopes. However, Ter 5 has chance to be detected by the 
VERITAS, and HESS telescopes close to ~1 TeV, i.e. at 
energies where the sensitivities of these telescopes are the 
largest. The eventual upper limits on the 7-ray flux from 
Ter 5 should constrain the acceleration efficiency of leptons 
inside this GC to the level of about one order of magnitude 
higher than in the case of 47 Tuc. Therefore, it seems that 
Ter 5 is not so promising GC from which GeV- TeV 7-ray 
emission might be observed in spite of the largest number 
of detected MSPs. Ter 5 is well situated at the large zenith 
angle for observations with the VERITAS telescopes and at 
low zenith angle for observations with the HESS telescopes 
what allows the coverage of a broadest possible energy range. 
Therefore, simultaneous observations at different zenith an- 
gles can provide useful information on possible variable 7- 
ray emission in the broad energy range. Note that discussed 
here model predicts steady 7-ray fluxes which should be eas- 
ily distinguished from possible variable 7-ray fluxes expected 
in case of their origin e.g. in the low mass X-ray binaries. 

5.3 M15 (NGC 7078) 

M15 is at the distance of 9.4 kpc and belongs to the class 
of the core collapsed globular clusters. M15 is very compact 
with the core radius of 0.07 arc min (R z = 0.19 pc) and a 
half mass radius of 1.06 arc min (Rh = 2.91 pc). The mass 
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Table 1. The upper limits on the parameter N p ■ r], describing the injection of leptons into the globular cluster, based on the sensitivities 
of different telescopes (G - GLAST; M - MAGIC; H - HESS; V - VERITAS). 



Model 


47 Tuc 


47 Tuc 


Ter 5 


Ter 5 


Ter 5 


M 15 


M 15 


M 15 


M 13 


M 13 


M 13 




(G) 


(H) 


(G) 


(M) 


(H, V) 


(G) 


(M) 


(H, V) 


(G) 


(M) 


(V) 


100 GeV, 2.1 


0.5 


0.01 


1.5 


0.8 


0.1 


0.6 


0.6 


0.05 


0.8 


0.5 


0.03 


100 GeV, 3.0 


0.15 


0.03 


1 


0.5 


0.3 


1.3 


0.4 


0.15 


0.5 


0.3 


0.1 


1 GeV, 2.1 


0.1 


0.03 


2 


2 


0.15 


0.5 


1.5 


0.1 


0.3 


0.8 


0.05 


1 GeV, 3.0 


0.15 


2 


0.6 


30 


6 


0.5 


20 


10 


0.5 


0.1 


3 


mono: 1 TeV 


1 


0.01 


1 


0.5 


0.02 


3 


0.5 


0.03 


2 


0.5 


0.03 


mono: 10 TeV 


2 


0.015 


3 


0.3 


0.05 


5 


0.8 


0.1 


5 


0.8 


0.15 



of the cluster is 5.4 x 10 5 Mq and the mass of the core inside 
10 4 AU is 3.4 x 10 3 M Q . It is supposed that the center of 
this globular cluster contains an ~ 2000 Mq mass black hole 
(e.g., Gerssen et al. 2000). M 15 contains 7 known millisec- 
ond pulsars. However, the existence of about 100 millisecond 
pulsars is expected (Tavani 1993). Moreover, two interest- 
ing binary systems have been observed in M15. One of them, 
belonging to the LMXB class, X2127+119, is observed in X- 
rays with the power ~ 10 36 erg/s (Beppo-Sax - Sidoli et 
al. 2000; Chandra - White et al. 2001). This system also 
shows huge outbursts with the peak power 6 x 10 38 erg/s 
(Samle 2001). The second one is an ultra-compact binary 
with on orbital period of ~ 17 min, (Dieball et al. 2005). 
The Whipple group searched for WIMP annihilation signa- 
tures which might surround the core of M15. The upper limit 
on the rate of the 7-ray photons, ~ 0.2 min" 1 can already 
weakly constrain some models of the WIMP mass and its 
distribution within M15 (Lebohec et al. 2003). 

We have calculated the 7-ray spectra expected from 
M15 in terms of the considered here injection scenarios of 
leptons. Due to comparable distances and masses, they are 
very similar to these ones obtained from Ter 5. The strongest 
constraints on the number of MSPs and the energy conver- 
sion efficiency (N p ■ 77) are also introduced by the telescopes 
able to observe at TeV 7-ray energies, provided that leptons 
are accelerated to such energies (see, Table 1). However, this 
GC is much better localized for the low zenith angle observa- 
tions with the MAGIC and VERITAS telescopes. Therefore, 
it is more interesting when one considers the detection of the 
7-ray signal in the case of the models which postulate accel- 
eration of leptons with relatively steep spectra. M15 can be 
also observed at large zenith angle with the HESS telescope 
providing additional constraints on the spectrum of injected 
leptons and the 7-ray variability in the broad range. The 
GLAST observations in the GeV energy range can also put 
interesting limits on energy conversion efficiency (see Table 
1) on the level below the present expectations from the the- 
oretical modeling of pulsars. For example, if 100 MSPs are 
present, then energy conversion efficiency from the pulsar 
wind to relativistic particles can be limited by the GLAST 
sensitivity to 77 < 0.005 - 0.013 (Table 1). Note that the ob- 
servations of the Crab Nebula limits 77 to ~ 0.1. Therefore, 
this parameter can be more significantly constrained by the 
observations of MSPs in globular clusters. 

Since, in principle, compact binary systems might also 



contribute to the high energy 7-ray emission from GCs, si- 
multaneous observations in broad energy range will allow 
extraction of the non-variable part of 7-ray signal expected 
in terms of the present model. 

5.4 M13 (NGC 6295): 

M13 is at the distance of ~7 kpc and belongs to the class 
of normal globular clusters. It has a core radius of 0.78 arc 
min (R c = 1.6 pc) and a half mass radius of 1.49 arc min 
(Rh = 3.05 pc). The total mass, 6 x 10 5 Mq, is comparable 
to M15. Already 5 millisecond pulsars have been observed 
within this cluster (periods 2-10 ms) and the existence of 102 
millisecond pulsars is predicted by Tavani (1993). The upper 
limit on the 7-ray emission, 1.08 x 10 -11 ph. cm -2 s _1 above 
500 GeV, has been reported by the Whipple collaboration 
(Hall et al. 2003). 

The 7-ray fluxes from M13 calculated in terms of the 
present model are shown in Fig. [5] These fluxes are a fac- 
tor of ~2 larger than in the case of M15, allowing to put 
stronger constraints on the number of MSPs and their en- 
ergy conversion efficiency (see, Table 1). M13 can be ob- 
served at small zenith angles by the MAGIC and VERITAS 
telescopes. Hence, it is good candidate for detection in the 
sub- TeV energies in case of injection of leptons with rela- 
tively steep spectra or with energies lower than ~ 1 TeV. 
The GLAST and MAGIC telescopes should be able to con- 
strain N p ■ r) on the level of <0. 1-1.0 (see Table 1). The 
constraints on N p ■ 77 introduced by the VERITAS observa- 
tion can be an order of magnitude lower. Therefore, for 100 
MSPs present in this GC, the energy conversion efficiency 
can be limited to ~ lO"' 3 " 4 '. Note, that the present upper 
limit put by the Whipple collaboration allows to constrain 
N p • 77 below ~1.6, which for supposed number of 100 MSPs 
inside M13 allows to put the limit on 77 0.016, which is 
about two orders of magnitude higher. 

5.5 Leptons injected directly from the inner 
pulsar's magnetospheres 

In order to check whether the population of millisecond pul- 
sars in specific GCs can produce 7-ray fluxes potentially ob- 
served by the mentioned above telescopes, we also perform 
calculations of the 7-ray spectra assuming mono-energetic 
injection of leptons from the pulsar magnetospheres. These 
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Figure 6. The differential gamma-ray spectra (multiplied by the 
energy square, SED) expected from 47 Tuc (NGC 104) for dif- 
ferent parameters of the injected spectrum of leptons: power-law 
spectrum with index a = 2 (solid curve), and 3 (dashed). The 
spectrum has low energy cut-off at E m i n = 1 GeV (thin curves) 
or 100 GeV (thick). 7-ray spectra produced by leptons with the 
upper energy cut-off at -B ma x = 30 TeV are shown on the up- 
per figure and at i? max = 3 TeV on the lower figure. The spec- 
tra are shown for the energy conversion efficiency from the pul- 
sar winds to leptons equal to t) = 0.01 and assuming the pres- 
ence of JV P = 100 millisecond pulsars inside the globular cluster 
(Np ■ r) = 1), with the parameters: the rotational periods equal 
to 4 ms and the surface magnetic fields equal to 10 9 G. These 
pulsar parameters determine the power of the pulsar winds de- 
scribed by Eq. [16] The available EGRET (E) upper limit for 47 
Tuc (Michelson et al. 1994) is marked by a double dashed line. 
The sensitivities of different 7-ray telescopes (5a- detection within 
50 hrs, Lorenz 2001) are marked by a double dotted curve (C 
- CANGAROO, H - HESS, and V - VERITAS), a double dot- 
dashed curve (M - MAGIC, and expected similar sensitivity of 
the HESS II in this energy range), and a double solid curve (G - 
GLAST). 



leptons can be accelerated inside the pulsar light cylinder 
radius as expected in the models for pulsed 7-ray emission, 
e.g, the extended polar gap model recently considered by 
Harding, Usov & Muslinov (2005). The power in relativistic 
leptons, L c ± , which leave the inner magnetosphere through 
the light cylinder radius, is determined by the rotational en- 
ergy loss rate of the millisecond pulsars, L p . Applying the 
recent model for acceleration of leptons inside the pulsar 
magnetopsheres, Muslinov & Harding (1997) estimates this 



O.ILp, where k ; 



0.15. 



power on, L e ± = 0.75k(1 — k)L p 
The Lorentz factors of leptons escaping from the magnetop- 
sheres of MSPs are expected to be limited by the curvature 
energy losses to values of the order of a few to several TeV 
(e.g. Bulik, Rudak & Dysk 2000; Luo, Shibata & Melrose 
2000; Harding & Muslinov 2002). In our calculations we ap- 
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Figure 7. As in Fig. 6 but for Ter 5. 
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Figure 8. As in Fig. 6 but for M 15. 



ply the characteristic values for leptons injected into the 
GCs of the order of 1 - 10 TeV. 

For easier comparison with the above calculations, 
which assume acceleration of leptons at the pulsar wind 
shock with the power-law spectrum, we also assume the pres- 
ence of 100 MSPs in specific GC and the efficiency of lepton 
acceleration equal to 77 = 0.01. The results of calculations 
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Figure 9. As in Fig. 6 but for M 13. 
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Figure 10. The differential gamma-ray spectra (multiplied by 
the energy square, SED) expected from 47 Tuc (solid curve), Ter 
5 (dashed), M15 (dot-dashed), and M13 (dotted) in the case of 
mono-energetic leptons injected from the pulsar inner magneto- 
sphere with energies 10 TeV (a) and 1 TeV (b). The 7-ray spec- 
tra are shown for the energy conversion efficiency from the pulsar 
winds to leptons equal to t) = 0.01 assuming the presence of 100 
MSP inside the GC. The sensitivities of different 7-ray telescopes 
and the EGRET upper limit are marked by letters as in Fig. [6] 



for 47 Tuc, Ter 5, M15 , and M13, in the case of mono- 
energetic injection of leptons with energies 1 TeV and 10 
TeV are shown in Fig. 1101 The 7-ray spectra expected from 
different GCs have similar shape but differ essentially in the 
intensity. Note interesting double-peak structure evident in 
these spectra in the case of 10 TeV leptons. They appear due 
to scattering of the MBR in the T regime (lower bump) and 
the scattering of stellar photons in the KN regime (narrower 
high energy peak). As in the case of the power-law spec- 
tra, the best candidates among GCs for positive detection 
with the Cherenkov telescopes are 47 Tuc (on the southern 
hemisphere) and M13 (on the northern hemisphere). Detec- 
tion of the GeV 7-ray flux by the GLAST telescope in the 
case of mono-energetic injection of leptons seems to be less 
likely. However, if the acceleration efficiency of leptons is 
one order of magnitude larger than considered in this sam- 
ple calculations, e.g. r\ ~ 0.1 as predicted by the inner pulsar 
acceleration models (Harding, Usov & Muslinov 2005), then 
also the GLAST telescope should detect the GeV 7-rays 
at least from 47 Tuc and M13. In case of non-detection of 
any signal in the 7-ray energy range, the present and fu- 
ture telescopes can put strong limits on N p • 77 (see Table 1). 
These limits allow to constrain energy conversion efficiency 
from the pulsar to relativistic leptons below r\ ~ 0.01 — 0.05 
(based on the GLAST observations in GeV energy range), 
»7~(3- 
and rj n 
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I x IQ~ A (the MAGIC observations at ~ 100 GeV), 
(the HESS and VERITAS observations at 



"(3+4) 



~ 1 TeV), provided that 100 MSPs are present in a specific 
GC (see Table 1). 



6 DISCUSSION AND CONCLUSIONS 

Globular clusters are promising sites for constraining the 
acceleration processes of leptons in the pulsar magnetop- 
sheres and/or at the shock waves in the relativistic pulsar 
winds. Expected large injection rates of leptons by the whole 
population of millisecond pulsars inside GC have to inter- 
act with strong thermal radiation field, created by the huge 
concentration of stars, and with the microwave background 
radiation. Therefore, it is likely that GCs should be accom- 
panied by a well localized, concentric, non-variable 7-ray 
sources. Here, we consider two scenarios for injection of lep- 
tons and calculate expected 7-ray fluxes in the GeV- TeV 
energy ranges from a typical globular cluster. A few specific 
sources (47 Tuc, Ter 5, M15, M13), in which large population 
of millisecond pulsars is expected, are considered in detail. It 
is concluded that all these objects should be detected either 
by the GLAST (in the GeV energies) or Cherenkov tele- 
scopes (in the TeV energies) depending on the shape of the 
spectrum of injected leptons provided that: (1) leptons are 
accelerated up to ~1 TeV; (2) specific GCs contain at least 
several MSPs; and (3) the power conversion efficiency from 
the pulsar winds to relativistic leptons is comparable to that 
one observed in the case of the Crab Nebula (rj ~ 10%). It 
is expected, that in the case of acceleration of leptons in the 
inner pulsar magnetosphere's, r\ is also of similar order, i.e. 
~ 3% (Harding, Usov & Muslinov 2002). 

In the case of lack of positive detection of 7-rays from 
GCs, our calculations allow us to put strong constraints on 
the product of the power conversion efficiency times num- 
ber of MSPs in specific GC (N p -77). In Table. 1 we show the 
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upper limits on iV p • rj for specific GCs which may be intro- 
duced based on the sensitivities of the present Cherenkov 
telescopes and the future GLAST telescope for a few dif- 
ferent models describing injection spectrum of leptons. If 
the acceleration of leptons occurs above ~ 1 TeV, then the 
most restrictive limits obtained by the HESS and VERITAS 
type Cherenkov telescope arrays on the power conversion ef- 
ficiency of leptons are of the order of r\ ~ lO - ' 3 " 4 ^, provided 
that N p ~ 100 MSPs are present in specific GC as expected 
in some estimations (e.g. Tavani 1993). The GLAST and 
MAGIC telescopes can provide also an order of magnitude 
more restrictive limits on r) than estimated from the observa- 
tions of the Crab Nebula or expected in the case of isolated 
MSPs, provided that leptons are accelerated by N p ~ 100 
MSPs to energies at least ~1 TeV (Table. 1). These limits 
put strong constraints on the high energy processes occur- 
ring in the inner pulsar magnetopsheres and/or the acceler- 
ation mechanisms in the pulsar winds and the pulsar wind 
shocks. 

Note that considered here scenario differs from those 
ones which postulate production of high energy 7-rays in the 
radiation processes occurring inside the inner pulsar magne- 
topsheres as proposed by e.g. Bulik, Rudak & Dysk (2000), 
Luo, Shibata & Melrose (2000), and Harding, Usov & Musli- 
nov (2005) or inside the interacting binary millisecond pul- 
sars (Tavani 1993). In one of discussed here scenarios, in 
contrast to earlier model considered by Tavani (1993), there 
is assumption that acceleration of leptons occurs at the part 
of the pulsar wind which is not limited by the wind of the 
companion star. The second considered scenario with the 
mono-energetic injection of leptons, in contrast to the inner 
magnetosphere 7-ray production models mentioned above, 
takes into account the radiation processes from relativistic 
leptons which escaped from the inner pulsar magnetospheres 
and propagate in the whole volume of the GC. 7-rays pro- 
duced in terms of those earlier models can additionally con- 
tribute to the expected 7-ray flux mainly in the GeV energy 
range. Note that also other mechanisms of 7-ray produc- 
tion can operate inside the GCs. For example, it can occur 
that GCs contain microquasar type X-ray binaries (recently 
established as TeV 7-ray sources), intermediate mass black 
holes, or related to them concentrations of dark matter parti- 
cles. In fact, powerful X-ray binaries have been also observed 
inside some GCs, e.g. mentioned above X2127+119 inside 
M15 and some evidences of the presence of the intermediate 
mass black holes inside GCs have been also reported, e.g. in 
M15 (Gerssen et al. 2002). Therefore, the constraints on the 
number of MSPs and the power conversion efficiency from 
the pulsar winds to energetic leptons derived in this paper 
should be considered as the upper limits. 
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